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Abstract 


Charmonium production from the hadron gas formed in ultra-relativistic 
heavy-ion collisions is studied. Using the J/if- hadron absorption cross section 
determined from the nonperturbative quark-exchange model, which has a 
peak value similar to that used in the comover model for J/if suppression 
and a thermally averaged value consistent from that extracted from J/if data 
in heavy-ion collisions, we find that J/ip production from the hadron gas 
is negligible in heavy-ion collisions at RHIC energies but is important at 
LHC energies as a result of the large number of charm mesons produced at 
higher energy collisions. The number of J/b produced from these secondary 
collisions at LHC may be comparable to that of primary J/i/j's, which are 
expected to be dissociated in the quark-gluon plasma created in the collisions, 
leading thus to a possible absence of J/^ suppression. Similar results are 
obtained for if' production in ultra-relativistic heavy-ion collisions. 
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One of the signals proposed for identifying the existence of a quark-gluon plasma in 
ultra-relativistic heavy-ion collisions is the suppression of J/'ip production compared to that 
expected from the superposition of nucleon-nucleon collisions [1]]. According to Matsui and 
Satz, if a quark-gluon plasma is created in heavy-ion collisions, J/ip's produced from initial 
nucleon-nucleon interactions will dissociate as a result of Debye screening and vanishing 
string tension between c and c. Recent experimental results from heavy-ion collisions at 
CERN SPS energies |2|-£j| have indeed shown a reduction of J /?/; production. Part of the 
effects can be attributed to absorption by nucleons as it is more likely that J/ip is first pro¬ 
duced as a pre-resonance in a color octet state and thus has a larger interaction cross section 
with nucleon |5|-||]. Such an effect is needed to account for the observed J/ip suppression 
in proton-nucleus collisions U, where one does not expect the formation of a quark-gluon 
plasma. Also, J /ip absorption by co movers, which are most pions produced in the collisions, 
has been suggested [TTi| [T3||. In collisions of light ions such as S+U at 200 GeV/nucleon, 
these absorption mechanisms are sufficient to explain the experimental data. However, for 
central collisions of heavy nuclei such as Pb+Pb at 160 GeV/nucleon, the anomalous large 
suppression of J/ip observed in the experiments has led to the suggestions that a quark-gluon 
plasma is formed in these collisions, and the J/ip yield is reduced due to dissociation either 
by Debye screening P~1|-[I6[] or collisions with gluons |T7|. If the quark-gluon plasma has 
already been formed in heavy-ion collisions at SPS energies, it is almost certain that it will 
also be formed in heavy-ion collisions at the Relativistic Heavy-Ion Collider (RHIC) at the 
Brookhaven National Laboratory and the Large Hadron Collider at CERN (LHC), where 
energies are much higher. Therefore, -J/v' suppression at RHIC and LHC is considered as 
one of the most prominent signals for the quark-gluon plasma. 

But will J/pj production at RHIC and LHC be really reduced? This depends on whether 
J/ip can be regenerated from the hadronic matter after the phase transition of the quark- 
gluon plasma. Such an effect has been shown to be important for hadrons made of strange 
quarks |TS|. In hadronic matter, J/ip can be produced from the interactions of charm 
mesons in reactions such as DD*,D*D,D*D* —> J/xpn. These reactions are apparently 
unimportant at SPS energies but may become significant at RHIC and LHC since charm 
production in nucleon-nucleon interaction increases with center-of-mass energy while the 
J/'ip to cc ratio remains essentially at a constant value of ~ 2.5 x 10 -2 
the reaction is exothermic as mp + m ^ ~ 3.73 GeV/c 2 and mo* + mg* 
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Furthermore, 
3.87 GeV/c 2 

while mj/^ + rn^ ~ 3.15 GeV/c 2 , it is thus more likely that J/^’s produced from the 
hadronic matter will survive and be detected in experiments. If this is the case, then the 
final number of J/'ip may be comparable to or even larger than that of primary J/'ip’’ s which 
are dissociated in the quark-gluon plasma, leading instead to an absence of J/'ip suppression 
or even an enhanced J/ip production. In the following, we shall make an estimate of J/'ip 
production from a pion gas formed in heavy-ion collisions at RHIC and LHC energies. 

To model heavy-ion collisions at such high energies, we use the results from the HIJING 
calculation |20|, which takes into account parton production from semihard scatterings. 
It shows that at an initial proper time r 0 a thermally equilibrated although chemically 
non-equilibrated quark-gluon plasma of temperature T 0 is formed. For Au+Au collisions, 
HIJING predicts that To ~ 0.7 and 0.5 frn/c, T 0 ~ 0.57 and 0.83 GeV at RHIC and LHC 
energies, respectively. The quark-gluon plasma then cools due to expansion and production 


of additional partons. It reaches the critical temperature T c ~ 


200 MeV at about r r 
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fm/c at RHIC and 6 fm/c at LHC, when the quark-gluon plasma starts to make a transition 
to a hadron gas, consisting mostly of pions. Since neither the partons have reached chemical 
equilibrium nor the order of phase transition is definitely known, the time for the system 
to remain at T c can not be properly determined. We shall assume that the proper time 
at which the quark-gluon plasma is completely converted to a pion gas is Th ~ 2r c and 
will study how the results depend on the value of Th- Since transverse expansion is not 
expected to be appreciable in this early phase, the initial volume of the hadronic matter is 
simply V h TTR^Th if one uses the boost-invariant model of Bjorken |H). In the above, i? 0 
is the radius of the colliding nuclei. Assuming that the hadron gas expands isentropically, 
its temperature then decreases according to the inverse of the cubic root of volume. To 
include the transverse expansion of the hadron gas, we introduce an acceleration a. Then 
its transverse radius increases with time according to R(r) = Ro + a{r — Th ) 2 /2 until the 
velocity reaches the velocity of light, when the transverse radius increases linearly with 
time. From entropy conservation, the temperature of the hadron gas then decreases as 
T{t) = (jh/ r) 1 / 3 (i? 0 / R{t)) 2 ^T c . Final hadron yields and spectra are determined at freeze 
out with a temperature Tf ~ 120 MeV. We note that for pions in equilibrium the density is 
approximately given by n n ~ 0.25(T/197MeV) 3 frn -3 . 

The time evolution of the temperature and volume of the hadron gas for Au+Au collisions 
at both RHIC and LHC is shown in Fig. |l] for a = 0.1 c 2 /fm. The transverse expansion 
velocity at freeze out is about 0.75c and 0.9c at RHIC and LHC, respectively. The initial 
number of charm quarks produced in these collisions can be estimated from their production 
cross section a^ 5 N in nucleon-nucleon interaction using [|22 

N £- B « ABa™/*™, (1) 

where cr^ 5 is the nucleus-nucleus (A+B) inelastic cross section. For central collisions, ctA b 
is approximately the nuclear geometrical cross section. The cross section <j^ N has been 
evaluated using the perturbative QCD [^0|, and it has a value 0.16 mb and 5.75 mb for 
nucleon-nucleon interaction at RHIC and LHC energies, respectively. For Au+Au collisions, 
the initial number of cc (i.e., DD) pairs, is thus about 3.8 at RHIC and 136 at LHC. To 
determine the initial J/Y number, we assume that the same cc to J/Y ratio, about 40, 
observed in nucleon-nucleon interaction at center-of-mass energy below 40 GeV |T9[ holds 


at RHIC and LHC energies. In this case, the primary J/ifj number is about 0.095 and 3.4 at 
RHIC and LHC, respectively. Since a quark-gluon plasma is almost certain to be created in 
heavy-ion collisions at RHIC and LHC, both J/^’s and D’s are expected to be dissociated 
due to Debye screening. As shown in Ref. |2(|], thermal charm production from the quark- 
gluon plasma is small, it is thus reasonable to assume that the c and c quark numbers remain 
constant during the quark-gluon phase. After phase transition to a hadron gas, these charm 
quarks are more likely to form D and D mesons and their excited states due to the large 
abundance of light quarks in the plasma. Conservation of charm then requires that the 
number of charm mesons at the beginning of the hadron phase is the same as the initial one. 

The J/ij} production cross section from DD interaction, can be related to the 

absorption cross section (J.j/^^dd via the detailed balance relation 

CDD^j/i/m — d(kjf^T T /k D £)Yaj^ 7r ^ D £ ) , ( 2 ) 

where the factor d is 3/4 for DD* and D*D annihilation, and 1/4 for D*D* annihilation; 
kj/^n and k D p are, respectively, the relative momenta of J/ijjit and DD. The magnitude 
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of J/ij) absorption cross section by pion, (Tj/^^od, is still under debate. In the comover 
model for J/ip suppression in heavy-ion collisions at SPS energies, it is taken to be about 
3 mb. Both perturbative QCD [6j and effective Lagrangian |^3[ calculations give a value 
which is about a factor of 10 smaller. On the other hand, studies including nonperturbative 
effects via quark-exchange model p4[ give a peak value of about 6 mb. As shown below, 


the thermally averaged effective cross section that takes into account the threshold effect is 
much smaller and is consistent with that extracted in Ref. 0 from J/-0 suppression data 
in heavy-ion collisions. In the following, we shall use the parameterized values from Ref. 


, i.e., 


0/#^dd(«) « (To (l - (y) 2 ) exp[—a(Vs - v/sl)]0(s - s 0 ), 


( 3 ) 


where s is the square of the J/ipi r center-of-mass energy and so is the threshold energy. For 
the final states DD* and DD*, the parameters are Co = 2.5 mb, s o = 15.05 GeV 2 , si = 17.6 
GeV 2 , and a = 11 GeV -1 ; while for the final state D*D*, they are ao = 3.4 mb, So = 16.2 
GeV 2 , si = 19.0 GeV 2 , and a = 11 GeV -1 . 

The quantity needed for estimating J/ijj production from DD interaction in a hadron 
gas is the thermal average of <Tdd^j/- 4 >k v , where v is the relative velocity between D and D, 

i.e., 


p OO 

<o D d^j/^v> = [4:(in D jT) 2 (m D jT) 2 K 2 (m Dl /T)]K 2 (rn D2 /T)]- 1 / dz 

J(m Dl +m D2 )/T 

■ [z 2 - ( m Dl /T + m D2 /T) 2 ][z 2 - ( m Dl /T - m D2 /T ) 2 ] 

• Ki(z)a D D^j/^n- ( 4 ) 

In the above, rrif 2l and mr> 2 are the masses of the two charm mesons; K\ and K 2 are, 
respectively, the modihed Bessel function of the first and second kind. The above quantity 
has been evaluated separately for DD*, D I). and D*D* annihilation. In Fig. |j. we show 
the temperature dependence of < odd*v >, defined by 


< &DD*V > — < 0’DD*^J/ilm v > + < OD*D->J/ilm v >j (5) 

and < a D *D*v >=< >. Similarly, one evaluates the thermal average 

< cr j/^^ D £)V > for the production of DD*, D*D, and D*D*, and defines < ojj^v > 
as their sum. The temperature dependence of < uj/^v > is also shown in Fig. We see 
that < cr D *i)*v > is about a factor of three smaller than < cr D ^v > due to a smaller cross 
section. Furthermore, < cr D jjV > is larger than < (Jj/^v > at all temperatures. The value 
of latter ranges between 0.25 to 1.0 mb, which is much smaller than the cross section itself 
as a result of threshold effects so only pions with sufficient energy can destroy a J/'ijj. 

The rate of J/-0 production from the hadron gas is given by 

dR 

—j— =< <Jdd*v > nono-+ < cr D *£>* > n 2 D *— < oj/^v > rij/^n n , (6) 

where rip and ud* are the densities of D and D*, respectively, and their relative value is 
determined by assuming that they are in chemical equilibrium, i.e., 
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Up _ 1 + (1/3) exp((?7iD* - m D )/T) 
n D * 1 + 3exp(-(m£,* - m D )/T) 


( 7 ) 


In Eq. (S|), we have used the fact that the densities of charm and anti-charm mesons are the 
same. The density of J/ip is denoted by nj /</,. 

The total number of produced J/ijj is obtained by multiplying the above rate by the 
volume of the hadron gas and then integrating over time. We find that at RHIC the number 
of J/-0 produced from DD annihilation is only 0.005 and is more than an order-of-magnitude 
smaller than the number of primary J/i]). This is different at LHC as a result of the much 
larger number of DD present in the hadron gas. In Fig. |3j, we show by the long-dashed 
curve the J/ip number calculated with a transverse acceleration a = 0.1 c 2 /fm as a function 
of time in Au+Au collisions at LHC. It is seen that the J/-0 number increases with time 
and at freeze out is comparable to the number of primary J/-0, which is indicated by the 
left arrow. We note that calculations carried out for heavy-ion collisions at SPS energies 
also show a negligible effect on J/-0 production from DD annihilation in the hadron gas. 

For larger values of a, the hadron gas expands faster so its temperature drops more 
appreciably, leading to a reduction of the number of produced J/if; as shown in Fig. |3] for 
a = 0.15 and 0.20 c 2 /fm. However, the transverse expansion velocity in these cases reaches 
the velocity of light after a few fm/c, which seems unreasonable. On the other hand, with 
smaller values of a, the final J/^ number is increased as shown in Fig. |3] for a = 0.05 c 2 /fm. 

If pions are out of chemical equilibrium [^], J/-0 absorption may be enhanced. Using 
a pion chemical potential of 130 MeV, the final J/ij) number is reduced only slightly and 
is still comparable to the primary one. The lifetime of the phase transition also affects the 
result. If it is doubled, then the number of produced J/ij) is reduced to about 2.25, while 
it is increased to 6.15 if the phase transition is instantaneous. Reducing the value of the 
temperature at which the phase transition occurs does not change much the number of J/ip 
produced from the hadron gas. Using T c = 150 MeV, the J/-0 number is reduced only 
slightly to 3.53 due to the competing effects of increased < <7dd* > and shortened lifetime 
of the hadron gas with decreasing temperature. 

The dependence of J/i^ production on the cross section (Jj/^^pp is more significant. 
If the cross section is an order of magnitude smaller than the one used here as predicted 
by the perturbative QCD 0 and the effective Lagrangian |23|] calculation, then the number 
of J/ij> produced from the hadron gas at LHC is only 0.41 and is much less than that of 
primary ones. It is thus important to have a better determination of this cross section from 
both experiments and theoretical models. 

The above analysis can be generalized to -0' production in ultra-relativistic heavy-ion 
collisions. The primary 0' number is about a factor of 5 smaller than the 0 number, i.e., 
about 0.019 and 0.682 at RHIC and LHC, respectively. Since the radius of 0' is about twice 
that of 0, its absorption cross section by a pion is about four times larger. Using this cross 
section and the 0' mass, we obtain only 0.002 0' from DD annihilation in hadron gas at 
RHIC, again negligible as in the case of J/if). For heavy-ion collisions at LHC, the number 
of 0' produced from the hadron gas is about 1.15 and is more than the number of primary 
0'. Therefore, there will not be a suppression of 0' production at LHC either. 

To summarize, we have estimated the number of J/0 produced from the reaction DD —> 
J/07T in the hadron gas formed in heavy-ion collisions at RHIC and LHC using the initial 
conditions determined from the HIJING parton model. Although this is a negligible effect 
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at RHIC, it may become important at LHC as a result of the appreciable number of D and 
D mesons in the hadron gas. We have found that the number of J/ip produced from this 
reaction at LHC is comparable to that produced from initial primary collisions, which are 
either absorbed by nucleons or dissociated in the quark-gluon plasma formed in the collisions. 
We thus conclude that in heavy-ion collisions at LHC one may not see a suppression of J/'ip 
production. A similar estimate has been made for ip' production, where the number of 
ip 1 produced from the hadron gas may be even larger than the primary ones. To use the 
J/ip and ip' yields as signals for the quark-gluon plasma at LHC therefore requires a good 
understanding of their production from the hadron gas. 
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T-xjj (fm/c) 


FIG. 3. Time evolution of the abundance of J/ij) and D from Au+Au collisions at LHC for 
different values of transverse acceleration a (c 2 /fm). The initial J/i^ number is denoted by the left 
arrow. 
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